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A hypothetical time-variation of the gravitational constant G would cause neutron star matter 
to depart from beta equilibrium, due to the changing hydrostatic equilibrium. This forces non- 
equilibrium beta processes to occur, which release energy that is invested partly in neutrino emission 
and partly in heating the stellar interior. Eventually, the star arrives at a stationary state in which 
the temperature remains nearly constant, as the forcing through the change of G is balanced by the 
ongoing reactions. Comparing the surface temperature of the nearest millisecond pulsar, PSR J0437- 
4715, inferred from ultraviolet observations, with our predicted stationary temperature, we estimate 
two upper limits for this variation: (1) \G/G\ < 2 x 10~ 10 yr -1 , if we allow direct Urea reactions 
operating in the neutron star core, and (2) \G/G\ < 4 x 10~ 12 yr" 1 , considering only modified Urea 
reactions. Both results are competitive with those obtained by other methods, with (2) being among 
the most restrictive. 

PACS numbers: 91.10.Op, 06.20.Jr, 97.60.Jd 



Since Dirac |l( suggested that the gravitational force 
may be weakening with the expansion of the universe, the 
question about whether or not the fundamental constants 
of nature vary in time has been of interest. General rela- 
tivity assumes a strictly constant gravitational coupling 
parameter G, but alternatives such as the Brans-Dicke 
scalar-tensor theory predict a variable G. 

There have been many experiments attempting to test 
the constancy of G, obtaining upper limits for its varia- 
tion, most usefully expressed as |G/G|. The experiments 
can be separated into three classes. In the first class 
are those that measure the variation of G from the early 
Universe to the present time, over timescales > 10 10 yr, 
such as Big-Bang nucleosynthesis and anisotropics in 
the cosmic microwave background radiation [4|, yield- 
ing upper limits on the long-term averaged variation, 
G/G| < 1CU 13 yr -1 . The second class contains measure- 
ments on long timescales, 10 9 ~ 10 yr, but without reach- 
ing into the very early Universe. In this category, one 
finds constraints based on paleontology []| and on stellar 
astrophysics, the latter including helioseismology j6j, bi- 
nary pulsar masses Q , and globular clusters || . Inferred 
constraints are in the range G/G < 10~( n ~ 12 ) yr -1 . 
Experiments in the third class measure the change of G 
directly over short, "human" timescales, ~ 10 yr. Orbits 
of planets [|J , rotation of isolated pulsars [h| i orbits of 
binary pulsars 01 > an d oscillations of white dwarfs 
are some examples of this class, with inferred upper lim- 
its |G/G| < l(H 10 - 12 ) yr" 1 . Even though results from 
the first category are nominally the most restrictive on 
a long-term variation of G, they depend on the assumed 
form of the variation of G near the Big Bang. Thus, it is 
still useful to consider measurements of the second and 
third categories, which could detect variations of G in 
more recent times. 



In this paper, we present a new method for constrain- 
ing a possible time- variation of G, which we call gravito- 
chemical heating. A change in G induces a variation in 
the internal composition of a neutron star, causing dissi- 
pation and internal heating. The comparison of the pre- 
dicted surface temperatures with the only available ob- 
servation 0| sets constraints on the variation of G. The 
method is based on the results of Fernandez & Reiseneg- 
ger ^3 (see also Refs. [l^, [l^, and \vi\). who showed 
that internal heating arises as a consequence of the com- 
pression of the star due to the decrease of its rotation 
rate and centrifugal force (rotochemical heating). In the 
present work, compression (or expansion) of the star is 
caused by an increase (or decrease) in the gravitational 
force as G changes, with the governing equations having 
the same mathematical form. Since the chemical relax- 
ation timescale of a neutron star is ~ 10 s yr, our ex- 
periment is closest to those in the second class discussed 
above. If G changes over cosmological timescales (as in 
the proposals of Dirac and Brans-Dicke) G can be ap- 
proximated as a constant, which we will do in most of 
what follows. 

Theory.- The simplest neutron-star models contain 
only neutrons (n), protons (p), and low-mass leptons (I; 
including electrons [e] and muons [p] ) , which can convert 
into each other by beta reactions: direct Urea, 

n — ► p + I + 77, 
P + l — ► n + y, (1) 

possibly forbidden by momentum conservation, and mod- 
ified Urea, 

n + N — > p + N + I +77, 
p + l + N — > n + N + v, (2) 
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where an aditional nucleon iV must be introduced in or- 
der to conserve momentum. Neutrinos (u) and antineu- 
trinos (v) produced in these reactions escape, taking en- 
ergy with them and therefore cooling the star 0] ■ When 
the latter is in beta equilibrium (with the chemical po- 
tentials satisfying ji n = fx p + for both I — e,/i), the 
rate of direct and inverse reactions is the same, i.e., the 
composition does not change. 

A change in G produces a change in the density of the 
stellar matter. Since the chemical potentials depend on 
density, the system departs from the beta equilibrium 
state. This departure can be quantified by the chemical 
imbalances, 



Vnpl = Mn — — Mi • 



(3) 



It increases the reaction rates so as to approach a new 
equilibrium configuration. If G changes continuously, the 
star is always out of equilibrium, storing an excess of en- 
ergy that is dissipated as internal heating and enhanced 
neutrino emission. 

The formalism for calculating the evolution of the tem- 
perature and chemical imbalances is described in § 2 of 
Ref. 0| . Here we outline the fundamental equations and 
the modifications required in order to treat the gravito- 
chemical case. The evolution of the internal temperature 
as measured by a distant observer, T (we avoid the usual 
subscript "oo" in order to keep the equations simple), 
taken to be uniform inside the star, is given by the ther- 
mal balance equation, 



f — — (L H - L v - L 7 ) , 



(4) 



where C is the total heat capacity of the star, Lh is the 
total power released by the heating mechanism (which 
depends on the chemical imbalances and the tempera- 
ture), L v the total power emitted as neutrinos, and L 1 
the power released as thermal photons. The evolution of 
the chemical imbalances is given by 

Vnpl = - [AD,l(Vnpe,T) + A M ,l{Vnpe,T)] 

- [BD,i(Vnp^,T) + B M ,i{rinpn,T)] + C n piG[5) 

where the subscripts D and M refer to direct and modi- 
fied Urea reactions, respectively. The functions A and B 
quantify the effect of reactions towards restoring chemi- 
cal equilibrium, and thus have the same sign of r\ np i |l4|. 
The constant C np i quantifies the departure from equi- 
librium due to G, being positive and depending on the 
stellar model and equation of state. It can be further 
decomposed as G np i = {Z npi - Z np )I G ^ + Z np I GtP , where 
Ia.i = (dN^ q I 8G)a is the change of the equilibrium num- 
ber of each particle species, N^ q , due to the variation of 
G, at constant total baryon number A, while Zi are coef- 
ficients depending only on the stellar model (see Ref. 0] 
for a rigorous calculation of the latter). 



As in the rotochemical heating case, the main conse- 
quence of gravitochemical heating is that the star arrives 
at a stationary state, where the rate at which G modifies 
the equilibrium concentrations is the same as the rate at 
which reactions drive the system toward the new equilib- 
rium configuration, with heating and cooling balancing 
each other 15]. The properties of this stationary state 
can be obtained by the simultaneous solution of equations 
(@J and JHJ with T = fj np i = 0. The existence of the sta- 
tionary state makes it unnecessary to model the full evo- 
lution of the temperature and chemical imbalances of the 
star in order to calculate the final temperature, since the 
stationary state is independent of the initial conditions 
(see Ref. [14J for a detailed analysis of the rotochemical 
heating case). For a given value of |G|, it is thus possible 
to calculate the temperature of an old pulsar that has 
reached the stationary state, without knowing its exact 
age. 

When only modified Urea reactions operate and nu- 
cleon Cooper pairing effects are negligible, it is possible 
to solve analytically for the stationary values of the pho- 
ton luminosity and chemical imbalances rj^ pl , as a 

function of stellar model and current value of |G/G|. The 
reason for this is that the longer equilibration timescale 
given by the slower modified Urea reactions yields sta- 
tionary chemical imbalances satisfying r\ np i 3> kT. The 
term Lh — L v in the thermal balance equation can be 
written asymptotically as KL e rj^ pe + KLfiVnpni where 
Klj are positive constants that depend only on stellar 
mass and equation of state . Assuming a stellar model 
with mass 1.41 Miscalculated with the equation of state 
A18+<5u+UIX* [l!j, the photon luminosity in the sta- 
tionary state comes out to be 



l;'~5x io 



2^ 



\G/G\ 



8/7 



10 



-12 



yr 



erg s 



(6) 



Assuming isotropic blackbody radiation, i.e., L 1 — 
4naR 2 T^, with R ~ 14 km the radius of the star as mea- 
sured by a distant observer, and a the Stefan-Boltzmann 
constant, the surface temperature of the star in the sta- 
tionary state is 



T? ~ 7.7 x 10 4 



\G/G\ 



2/7 



10 -12 yr -l 



K. 



(7) 



The timescale for the system to reach the stationary 
state, obtained by the same procedure as in § 4.4 of 
Ref. 0, is 



Tst ~ 2 x 10 8 



\G/G\ 



-6/7 



1Q -12 yj-1 



yr. 



(8) 



Results.- In order to apply this formalism and con- 
strain the value of |G/G|, we need a neutron star that 
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(1) has a measured surface temperature (or at least a 
good enough upper limit on the latter), and (2) is con- 
fidently known to be older than the timescale to reach 
the stationary state for the intended upper limit on 
\G/G\. So far, the only object satisfying both condi- 
tions is the millisecond pulsar closest to the Solar Sys- 
tem, PSR J0437-4715 (hereafter "J0437"), whose surface 
temperature was inferred from an ultraviolet observation 
by Kargaltsev, Pavlov, and Romani 13] . Its spin-down 
age, r s( j — 5 x 10 9 yr (e.g., [Ill), and the cooling age of its 
white dwarf companion, twd — 2.5 — 5.3 x 10 9 yr [22|, are 
much longer than the time required to reach the steady 
state, as long as \G/G\ > 10~ 13 yr -1 . 

Consequently, we consider stellar models constructed 
from different equations of state, with masses satisfy- 
ing the constraint obtained for J0437 by van Straten et 
al. HJ, M psr = 1.58 ± 0.18A/©, and calculate the sta- 
tionary temperature for each, ignoring potential Cooper 
pairing effects 18]. In Figure 2] we compare the ob- 
servational constraints on the temperature of this pulsar 
[l3| together with the theoretical prediction for differ- 
ent equations of state (for further details see Table 1 of 
Ref. 0), assuming \G/G\ = 2 x 1CT 10 yr" 1 . When the 
stellar mass becomes large enough for the central pressure 
to cross the threshold for direct Urea reactions, T s drops 
abruptly, due to the faster relaxation towards chemical 
equilibrium. This occurs in two steps, as electron and 
muon direct Urea processes have different threshold den- 
sities (see, e.g., Ref. Q). The chosen value of G is such 
that the stationary temperatures of all stellar models lie 
(just) above the 90 % confidence contour, as seen in Fig- 
ure ^ an d therefore represents a rather safe, general up- 
per limit. 

Nevertheless, conventional neutron star cooling mod- 
els reproduce observed temperatures better when only 
modified Urea reactions are considered Restricting 
our sample to the equations of state which allow only 
modified Urea reactions in the mass range considered 
here, namely A18 + Sv, A18 + Sv + UIX, BPAL21, 
and BPAL31, we obtain a more restrictive upper limit, 
\G/G\ < 4 x 1CT 12 yr~\ illustrated in Figure El This 
value is among the most restrictive results from methods 
probing a similar timescale. 

So far, we have assumed a steady, monotonic time- 
variation of G. Can our method constrain short-timescale 
variations? We consider an oscillatory variation, G(t) = 
Go[l + Acos(uit)], in a star with only modified Urea re- 
actions and one leptonic species (electrons), in the inter- 
esting limit rj npe 3> kT. The equation for the chemical 
imbalance evolution is 

Vnpe = -K A ri 7 npe + C npe G. (9) 

In the regime of low forcing frequency (when the period 
is longer than the timescale to reach the steady state), 
this evolution is equivalent to the situation already con- 
sidered. In the opposite limit, when the oscillation is 
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FIG. 1: Comparison of the gravitochemical heating predic- 
tions to observations of PSR J0437-4715. The solid lines are 
the predicted stationary surface temperatures as functions of 
stellar radius, for different equations of state (A18 from [l9| 
and BPAL from pHl - ). constrained to the observationally al- 
lowed mass range for this pulsar |20|. Dashed lines correspond 
to the 68% and 90% confidence contours of the blackbody fit 
of Kargaltsev et al. |l3l| for the ultraviolet emission from this 
object. The value of \G/G\ = 2 x 10~ 10 yr -1 is chosen so 
that all stationary temperature curves lie above the observa- 
tional constraints. (BPAL32 and BPAL33 allow direct Urea 
reactions in the observed mass range of J0437.) 
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FIG. 2: Same as Figure Q but now the value of |G/G| = 
4 x 1CP 12 yr -1 is chosen so that only the stationary tem- 
perature curves with modified Urea reactions are above the 
observational constraints. 



faster than the time required to reach the steady state, 
we have 77 npo — C n p e AGo cos(wi). As mentioned above, 
the photon luminosity satisfies L 7 ~ KLeiJnpe- Assuming 
a detection limit similar to the luminosity observed in 
J0437, Lj ~ 10 29 erg s _1 , the heating can be detected 
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only if the oscillation timescale satisfies 



through its Regular Grant 1060644. 



1 {L 1 jK Le 



Cnpc|G 



> 3 x 10 s 



\G/G\ 



10 



-12 



yr" 



yr, 



(10) 



too long to be of interest. 

Conclusions. - We have shown that a long-term, mono- 
tonic variation of G alters the chemical equilibrium state 
of a neutron star, leading to internal heating. This ef- 
fect is an exact analog of that produced by spin-down 
compression [l4|. As in that case, the star reaches a sta- 
tionary state in which the temperature remains nearly 
constant. We have calculated this temperature as a func- 
tion of \G/G\ and stellar mass using realistic equations 
of state, and compared our results to the temperature 
inferred from ultraviolet observations of the millisecond 
pulsar J0437-4715 by Kargaltsev et al. I n the most 
general case, when direct Urea reactions are allowed to 
operate, we obtain an upper limit \G/G\ < 2 x 10~ 10 
yr -1 . Restricting the sample of equations of state to 
those that allow only modified Urea reactions, we obtain 
a much more restrictive upper limit, \G/G\ < 4 x 10~ 12 
yr , competitive with constraints obtained from other 
methods probing similar timescales. However, since the 
composition of matter above nuclear densities is uncer- 
tain and millisecond pulsars are generally expected to be 
more massive than classical pulsars, we cannot rule out 
the result for the direct Urea regime. Information about 
the prevalence of the direct Urea process and the mass 
threshold above which it is present may in the future be 
obtained from the study of the quiescent emission of soft 
X-ray transients (e. g., Ref. [18j). 

Further progress in our knowledge of neutron star mat- 
ter will allow this method to become more effective at 
constraining variations in G. The method can also be 
improved with an increased sample of objects with mea- 
sured thermal emission. In this respect, an interest- 
ing candidate is PSR J0108-1431, plausibly the nearest 
pulsar, with distance estimates ranging from 60 pc |2^| 
to 184 pc 24]. Due to its much slower rotation (pe- 
riod P = 0.808 s [U, compared to P = 0.00576 s for 
J0437 0|), the current time-derivative of its centrifu- 
gal force (oc P/P 3 ) is 680 times smaller than that of 
J0437. Thus, unless it was born rotating much faster 
(Pq < 0.02 s), rotochemical heating should be negligi- 
ble for this pulsar 0, and any detectable surface 
temperature would be a clear signature of a change in G. 

This work was supported by FONDECYT (Chile) 



[1] 
[2] 
[3] 

[4 

[5 

[6 

ir 

[8 

[9 

[10 
[11 

[12 

[13 

[14; 

[15 
[16 
[17 

[is; 

[19 
[20 

[21 
[22 
[23 

[24 



[25 



P. Dirac, Nature (London) 139, 323 (1937). 

C. Brans and R. Dicke, Phys. Rev. Lett. 124, 925 (1961). 

C. J. Copi, A. N. Davis, and L. M. Krauss, Phys. Rev. 
Lett. 92, 171301 (2004). 

R. Nagata, T. Chiba, and N. Sugiyama, Phys. Rev. D 
69, 083512 (2004). 

W. Eichendorf and M. Reinhardt, Mitteilungen der As- 
tronomischen Gesellschaft Hamburg, 42, 89 (1977). 

D. B. Guenther, L. M. Krauss, and P. Demarque, Astro- 
phys. J. 498, 871 (1998). 

S. E. Thorsett, Phys. Rev. Lett. 77, 1432 (1996). 

S. Degl'Innocenti et al. Astron. Astrophys. 312, 345 

(1995). 

J. G. Williams, X. X. Newhall, and J. O. Dickey, Phys. 
Rev. D 53, 6730 (1996). 

I. Goldman, Mon. Not. R. Astron. Soc. 244, 184 (1990). 
V. M. Kaspi, J. H. Taylor, and M. F. Ryba, Astrophys. 
J. 428, 713 (1994). 

O. G. Benvenuto, E. Garcfa-Berro, and J. Isern, Phys. 
Rev. D 69, 082002 (2004). 

O. Kargaltsev, G. G. Pavlov, and R. Romani, Astrophys. 
J. 602, 327 (2004). 

R. Fernandez and A. Reisenegger, Astrophys. J. 625, 291 
(2005). 

A. Reisenegger, Astrophys. J. 442, 749 (1995). 

A. Reisenegger, Astrophys. J. 485, 313 (1997). 

A. Reisenegger, P. Jofre, R. Fernandez, and E. Kantor, 

Astrophys. J. , in press (astro-ph/0606322 ). 

D. G. Yakovlev and C. J. Pethick, Ann. Rev. Astron. 

Astrophys. 42, 169 (2004). 

A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, 
Phys. Rev. C 58, 1804 (1998). 

W. van Straten, M. Bailes, M. Britton, S. R. Kulkarni, 
S. B. Anderson, R. N. Manchester, and J. Sarkissian, 
Nature (London) 412, 158 (2001). 

M. Prakash, T. L. Ainsworth, and J. M. Lattimer, Phys. 
Rev. Lett. 61, 2518 (1988). 

B. M. S. Hansen and E. S. Phinney, Mon. Not. R. Astron. 
Soc. 294, 569 (1998). 

T. M. Tauris, L. Nicastro, S. Johnston, R. N. Manchester, 
M. Bailes, A. G. Lyne, J. Glowacki, D. R. Lorimer, and 
N. D'Amico, Astrophys. J. 428, L53 (1994). 
Calculated online from the Cordes-Lazio NE2001 
Galactic free electron distribution model, 
http://rsd-www.nrl.navy.mii/7213/lazio/ne _model/l 
(J. M. Cordes and ~. T. W! Lazio, preprint 
astro-ph/0207156). 

D. Gonzalez, A. Reisenegger, and R. Fernandez, paper in 
preparation (2006). 



